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Motivated by the possibility that the ampUtude for neutrinoless double beta de- 
cay may be much smaller than the planned sensitivity of future experiments, we study 
ansatze for the neutrino mass matrix with M^e = 0. For the case in which CP is 
conserved, we consider two classes of real-valued mass matrices: "Class I" defined 
by |Me^| = |MeT-|, and "Class 11" defined by = |Mt-t-|. The important phe- 

nomenological distinction between the two is that Class I permits only "small" values 
of VeS up to ~ 0.03, while Class II admits "large" values of Ves up to its empirical 
upper limit of 0.22. Then we introduce CP-violating complex phases into the mass 
matrix. We show that it is possible to have tribimaximal mixing with Mf>e = and 
|M^T-| = |-/Vf^^| = |Mt-t-| if the Majorana phase angles are ib7r/4. Alternatively, for 
smaller values of |-/Vf^T-| = |M^^| = |Mt-t-| it is possible to obtain |V^3| ~ 0.2 and 
generate relatively large CP-violating amplitudes. To eliminate phase redundancy, we 
emphasize rephasing any mass matrix with M^e = into a standard form with two 
complex phases. The discussion alternates between analytical and numerical but re- 
mains purely phenomenological, without any attempt to derive mass matrices from a 
fundamental theory. 

1 Data and Conventions 

The present empirical knowledge of neutrino oscillations can be summarized qualitatively 
as follows jH E]. We observe a deficiency of electron neutrinos originating from the sun 
and attribute this to oscillations described roughly by a mixing angle ^soiar ~ 0.5 — 0.8 
and a mass-squared difference /S.ml^-y^^ ~ 10~^ — 10~^ eV^. We also observe a deficiency of 
muon neutrinos in the earth's atmosphere from incident cosmic rays and attribute this to 
oscillations described roughly by a mixing angle 6'atm ~ 0.6 — 1 and a mass-squared difference 
Am^^-jjj ~ 10~^ eV^. The commonly accepted theoretical interpretation of the data is that all 
three flavors of neutrinos - t'e, v^, and - participate in oscillations. In this work we base 
our quantitative empirical understanding of three-flavor neutrino oscillations on the analysis 
of Gonzalez- Garcia and Maltoni [2], who report the entries of the 3-by-3 neutrino mixing 
matrix V as having magnitudes^ 



^ There are varying degrees of confidence levels assigned to the different fits in the review. To impose as 
little theoretical prejudice as possible, we will always quote the 3cr bounds, which are the least restrictive. 
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where the bounds are correlated such that V is unitary. We also quote the recently updated 
report by Gonzalez-Garcia, Maltoni and Salvado [3] for the angle^ 

0.620 < ^1 < 0.934 , 0.00 < 62 < 0.218 , 0.550 < ^3 < 0.658 

and the mass-squared differences 



- m? = 7.59 ( ^^^gg ) x 10'^ eV^ , and 



h2.46 ± 0.37 X IQ-^ eV^ ("normal hierarchy") 
m^-m^- <^ -2.36 ± 0.37 x lO^^ eV^ ("inverted hierarchy") ' 

Although we know the mass-squared differences m'fj = mf — m^, we do not know the actual 
value of any of the mj. Thus to compare with oscillation data, we compute the ratio of 
mass-squared differences 

_ mg — _ J +25.5 to +41.0 ("normal" hierarchy) , 
^= ml -ml ^ 1-39.6 to -24.3 ("inverted" hierarchy) ' ^ ' ^ 

To study the mixing matrix, we use the standard angular parameterization |1] for unitary 
matrices given by 

V = ICVpumM , where K. = diag{e"'\ e'^\ e'^') , M = diag(e*^ e^", 1) (1.3) 

and 



^MNS = 

—C2C3 C2S3 
C1S3 + S1S2C3 C1C3 - S1S2S3 S1C2 I . (1.4) 

K^lSs — C1S2C3 S1C3 + C1S2S3 —CiC2j 

Here cj = cos 6'/, sj = sm9j and §2 = 826^^'^^, and we have chosen the sign conventions in 
^MNS to minimize the number of minus signs that appear. The angles in /C are unphysical 
and can be chosen arbitrarily. 




We will assume that the neutrinos are Majorana. In this case the Majorana phase ma- 
trix Ai is physically meaningful and contributes to the amplitude for neutrinoless double 
beta decay. In addition, the neutrino mass matrix M is symmetric. We will work in the 
basis for which the charged lepton mass matrix is diagonal with real positive entries, called 
the flavor basis. In this basis the neutrino mass matrix is 

M = V*DV^ (1.5) 

^The report quotes two sets of ranges for the angles, depending on uncertainties in the capture cross 
section of gaUium. The distinction between the two sets is a shght change in the range of 9^ and in the 
upper bound of 6*2, with 9i and the mass-squared sphttings unaffected. We take the least restrictive bounds 
whenever possible. Also, the notation in the reference is 6123 = Oi, di^ = 62 and 6*12 = 63. 
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where D = diag(mi, 1712, 1TL3). Here mj > denote the physical masses of the three neutrinos. 

There is some degree of rephasing freedom in the neutrino mass matrix M, and we will return 
to this point in a later section on CP violation. For now we simply wish to clarify a potential 
source of confusion for the case in which CP is conserved. If CP is conserved, then M can be 
taken as real, and we can without loss of generality set 6cp = 0. However, we cannot set p 
and a equal to zero, since the Majorana phase matrix Ai appears squared in the mass matrix 

M. With /C = J and 6cp = 0, we have M = /C*1/p*mns-^*^-^*'^pmns^* = '^pmns^ Vp^MNS> 
where we have defined the diagonal matrix 



D = j m2 with nil = f^i e and ni2 = ni2 e '^"^ . (1-6) 

The notation is such that mi^2 are real and positive while mi 2 are complex. (Also, with our 
phase conventions is always real and positive.) Thus the choice of or | for p and a 
generates non-removable minus signs associated with mi and ni2, which yield qualitatively 
different textures for the mass matrix M. In the CP-conserving case, it is convenient to 
separate these signs from V and instead associate them with the diagonal matrix D. 

2 Neutrinoless Double Beta Decay and Mee = 

As discussed at the end of the previous section, the choice of signs in mi 2 imply qualitatively 
different textures for the mass matrix. To motivate a particular choice, we recall the well- 
known fact that a direct way to measure one of the entries in M is in neutrinoless double 



beta decay, the amplitude of which is proportional to |Mee|. From Eqs. (1.3 1.5) we have 



\M^e\ = \cl (c^e^^^mi + s^e^^^ms) +s^m3| (2.1) 

where ai = — 2(p + 6cp) and 02 = — 2((j + Sqp)- Thus in general, \Mee\ depends on all three 
masses mi,m2 and m3, the two angles 62 and 63, and two phases ai and a2- 

A brief review of the current status of neutrinoless double beta decay was given recently 
by Bilenky [2] , which we now summariz The Heidelberg-Moscow and CUORICINO exper- 
iments imply the upper bounds \Mee\ < (0.3 — 1.2) eV and \Mee\ < (0.3 — 1.7) eV, respectively. 
The future experiments CUORE, EXO, GENIUS and MAJORANA plan to significantly im- 
prove the sensitivity to roughly |Mee| ~ (1 — 7) x 10~^ eV. 



These values should be understood in comparison to the mf- data above Eq. (1.2), which 



imply 



8.3 X 10-^eV < ^Jml^ < 9.1 x lO^^eV and 4.6 x lO'^eV < y'mli < 5.3 x lO^^eV 
for the normal hierarchy mi < m2 < m3. 



^An early review of neutrinoless double beta decay was given by Zel'dovich and Klhopov [B]. 
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If ms ^ rni^2, then rris ~ ^/fn^i ~ 5 x 10^^ eV, but since < 4.68 x 10^^ the large ms is sup- 
pressed by the small in |Mee|. Thus |Mee| is at most ~ 10^^ eV, which is an order of mag- 
nitude smaller than the planned sensitivity of future experiments. If rris > mi^2 but all three 



masses are still almost equal, then the term drops out and 

1 below Ec 
ai and 62 



"^1^2 1 C3 + 53 e 



2 p j(a2-ai)| 



With the bounds given below Eq. (1.1), this implies 0.26 < 
bound occurs for a2 



Mee\/"mi < 1, where the upper 



0. 



Thus for any normal hierarchy, \Mee\ tends to be smaller than the other entries in M. Using 
this as guidance, we suppose that Mee could be tiny and thereby set M^e = 0. In other words, 
throughout this paper we assume that the amplitude for neutrinoless double beta decay is 
zero, at least as a leading order approximation [3 El El ITU] . 

3 Tribimaximal Mixing with Afge = 

As has been noted independently by many authors [HI |12], the theoretical ansatz of "tribi- 
maximal mixing" defined as 



PMNS 



Vtr = 





1 


0\ 






1 


1 


1 


7?/ 



-0.82 0.58 
0.41 0.58 0.71 
0.41 0.58 -0.71, 



(3.1) 



is compatible with the empirical bounds given in |V^xp|- If neutrino oscillations conserve CP, 
then we can write the neutrino mass matrix in the flavor basis as M = Vrmns-DI^mns ' 
where D = diag(mi, m2, m^) and rhi = ±mi with uncorrelated signs. We can thereby define 
a "tribimaximal mass matrix" Mtb = Vtb D associated with the ansatz of tribimaximal 
mixing. Explicitly, this mass matrix read^ 



TB 



1 

3 



2 -1 -i\ n 1 V 

mi I 1/2 1/2 + ma I 11 



1/2, 




(3.2) 



For all values of rhi and ms, this matrix exhibits the symmetry Mg^ 
[T3t fT7] . We stress that although the condition V = Vtb necessarily implies Mg^ 



Mer and M^^ = M 



TT 

and 



M, 



Mt-t-, the converse is not true: M^^ 



and M^^ 



M-rr do not necessarily imply 



tribimaximal mixing. 



Two appealing examples of tribimaximal mass matrices with Mee = are obtained by choos- 
ing the values (mi,m2,m3) = (—1,2,9) and (—1,2,11), which give 




and M 



(-1,+2,11) 



TB 




(3.3) 



respectively. The first has R ~ 27, while the second has i? = 40 exactly, which correspond 
nearly to the lower and upper empirical bounds for R. 



^Since the mass matrix is symmetric, we display explicitly only its upper triangle. 
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In the mass matrix Mtb the sign flip [7711,1712) — ?► (— wii, —7712) effects the exchange 

(M^^, M^r) ^ -{M^r, M^^) . (3.4) 

This means that given one tribimaximal mass matrix, we can always find a second tribimax- 
imal mass matrix by interchanging the magnitudes of M^^ and M^t-. 



Thus from (3.3) we can immediately write the matrice^ 




M^^'~^'^^ =m^\^ 4 -5 j and M^^'~^'^^^ = { 5 -6 | (3.5) 
which also predict R ^ 27 and R = 40 respectively. 

A mass matrix that resembles the examples given above but with non-tribimaximal mix- 

J6l 

■ 0.84 0.54 0.02 
R^29 and |1/pmns| ~ | 0.33 0.54 0.77 | (3.6) 

0.43 0.65 0.63 

which was suggested in the context of a particular model [15] . Since M^^ 7^ Mrr the resulting 
mixing matrix is not tribimaximal, as can be seen from the nonzero Ves- On the other hand, 
both MnTB and Mtb share the property M^^ = Mg^. The matrix Mqtb also has M^e smaller 
than the other entries. 

In an attempt to systematically study this distinction, we consider the "Class I" ansatz 

/O 1 1 \ 

Ml = 771 A M^^ (Class I) (3.7) 

with = and M^^ = M^t-. To further classify deviations from tribimaximal mixing with 
Mee = 0, we also consider the "Class 11" ansat^ 

/O Me A 

Mii = mJ 5 M^rj (Class II) (3.8) 

with M^^ = Mrr- Since oscillation experiments cannot determine the overall scale of M, we 
from now on set mj, = 1 and treat the entries of M as dimensionless numbers. 

We emphasize to the reader that we make no attempt to derive these mass matrices from 
any theoretical model but instead study these matrices on purely phenomenological grounds. 





^We have used the rephasing freedom in M to move around the minus signs. See Section 4 
^Many authors have proposed parametrizations of deviations from tribimaximal mixing [14j . 
^In Class II, the value A/^^ = ^fi-r = 5 is merely a convenient normalization for comparing the empirically 
allowed mass matrices with those of Class I. 
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4 Rephasing the Mass Matrix 



Before proceeding to study the matrices Mi and Mn, we should comment on the significance 
of various signs that may appear in the mass matrix. Consider the most general 3-by-3 
complex symmetric matrix M with Mge = 0: 



'0 aef.e''^-^ aerc'"^-- 



M= a^f^e'^^^ a^re'^P^^ (4.1) 



where a^p and ipap are real numbers. Using the form M = V*DV'^ with V = /CVpmns-^ 
introduced in Section [1} we have M = /C*M/C*, where M = ^pmns-^*-^-^*^pmns' ^^us 
M = JCMK,. We are free to choose the phases in K, as we please, since they are unphysical. 
Choosing Hi = ^^p^,^- ^Pe^l, 1^2 = -^Vfif^ and K3 = -^ip^r gives 



M = \ a^^ a^.e*^ I (4.2) 



where ip = ^eT — Ve^l + \v^l^^~\VTT and rj = fiiT — -^'Pnfi — \fTT- We may thus dispense with the 
matrix M and consider only the matrix M. Henceforth when there is no risk of confusion we 
put any mass matrix M into the form of M and then drop the hat for notational convenience. 

For the case in which M is real, the phases reduce to the signs ±1. The above argument 
shows that any real-valued neutrino mass matrix with Mee = can be put into the form 

/O I Me J 

M= \M^^\ C\M^r\ I (4.3) 

where each of ( and (' can be either +1 or —1. The matrix M can be multiplied on both 
sides by the matrix Z = diag(l, 1, — 1), which transforms (CO ~^ {~Cy~C') and thereby 
leaves the product CC' unchanged. Since det{ZMZ) = det M, all observables based on the 




M in (4.3) are invariant under M — )■ ZMZ and therefore depend only on sgn(MerM^T-), not 
on ( = sgn(Mer) and (' = sgn(M^T-) individually. 

If we allow Mf^T- to range over all real numbers, then in both Classes I and II we can take all 
other entries in M to be strictly non-negative. Given this choice, it will turn out furthermore 
that only M^,- < can fit data. This can be seen from the form of Mtribi with 1713 ^ mi^2- 

To summarize, we will first study real- valued mass matrices of Classes I and II given in 
Eqs. (3.7) and (3.8) with M^f, = 0, M^,- < and all other entries positive. 



5 Analytic Preliminaries 

If CP is conserved in the neutrino sector, there are 6 potential observables in neutrino phe- 
nomenology: 3 angles 6i and 3 masses mj. Accordingly, a general 3-by-3 real symmetric 
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matrix has 6 independent parameters and thereby makes no predictions. By fixing Mee = 
we impose a constraint and thus fix one of the parameters [15], namely the angle 62 (and 
thus Ves), according to the relation 



tan^e2 = -('^^^^±^) . (5.1) 

Empirically we know that 62 < 0.22, so tan^6'2 ^ 1. This tells us that we cannot have 
ms ^ mi ~ m2, thus forbidding the inverted hierarchjipl for the ansatz Mee = 



As a limiting case, for ^2 — ^ we predict mi — )■ m2 tan^ ^3 and thus fix all three neu- 
trino masses. Since 0.36 < tan^^g < 0.60 and 6.9 x lO^^eV^ <'m\-m\< 8.2 x 10"^ eV^ 
we have ^ 

lim mi = \ — \- = 6.4 x 10"^ eV to 1.1 x 10"^ eV 

62^0 V cot^ - 1 



and 



lim mg = 4.6 x 10"^ eV to 5.4 x 10"^ eV . 

6*2^0 



Note that these ranges are rather narrow: mi can be only as large as 11/6.4 ~ 1.7 of its 
minimum value, and can be only as large as 54/46 ~ 1.2 of its minimum value. 

After fixing Mee = 0, the next step is to specialize either to Class I by imposing Mg^ = Mer 
or to Class II by imposing M^^ = M^^. Either choice will fix the angle 61 in terms of the 
other parameters, thus reducing the number of free parameters to four: the three masses m^ 
and the angle 6*3. 

If we were to impose the condition Me^ = Mer (Class I), then this would fix 

1 - XS2 , c|mi + s|m2 - ma 

tant^i = where x = — ^— — . 

I + XS2 C3S3[m2 - mi) 

Instead, if we were to impose the condition M^^ = M^r (Class II), then this would fix]^ 
tan 6*1 = —y + a/?/^ + 1 where 

_ (mi -m2)s2sin(203) _ 2 22 1 _ 2 22 

y = z r — With p = — SoCo and g = — SnSo . 

021713 — p mi — q 1712 

In the limit 62 — )■ 0, the two conditions become equivalent and imply 61 — )■ tt/A. Therefore 
the limiting case ^2 ~^ of our matrices of Classes I and II corresponds to the /ir-symmetric 



*In more detail, the inverted case TO3 < mi < m2 would require roughly | ^ tan^ ^'3 ± 1| < 6 x 10~^, where 
the ± is fixed according to rhirh2 > (plus sign) or mim2 < (minus sign). For the + case, this inequality 
is clearly impossible to satisfy since 1 is larger than 10^^. For the minus case, the ratio 1112 /mi would have 
to be of order 1 but fine-tuned to two decimal places. We will not consider this particular case and thereby 
specialize to mi < m2 < m^. 

^This condition results in a quadratic equation for tan 61 whose two roots are tan 9i = —y ± -y/y^ + 1 '. 
We choose the + sign to keep tan 6*1 positive. 
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ansatz 



M, 



fiT—sym 




(5.2) 



which along with the possibihty of = 7r/4 was studied by many authors |17j . 



At this point we should comment on fir symmetry in the neutrino mass matrix. Since 
the mass of the r is an order of magnitude larger than the mass of the muon, the effective 
Lagrangian at energy scales below already exhibits deviations from any underlying fir 
symmetry that may exist at high energy. At the energy scale of neutrino masses m^, ^ m^, 
any high-energy /ir symmetry should be badly broken and thus corrections to M are to be 



expected in general. Thus the /zr-symmetric texture of (5.2) should be thought of at most 
as a useful starting point for a phenomenological analysis. 



It is also worth remarking that tribimaximal mixing implies M, 



and M^^ 



but the converse is not true. For the /ir-symmetric ansatz (5.2), the third eigenvector is 
Va3 oc (0,1,-1), exactly as for tribimaximal mixing, irrespective of the values of a and b. 
However, the other two eigenvectors are only proportional to (—2, 1, 1) and (1, 1, 1) for the 
particular case b = 1 — a. For example, a = 5 and 6 = 1 — 5 = —4 reproduce the matrix 
of Eq. (3.3). On the other hand, changing b to —4.2 gives 



M 




R^36 and IK 



PMNS 



0.80 0.60 
0.43 0.56 0.71 
0.43 0.56 0.71 



which deviates from tribimaximal mixing in the first two columns of V. Attempting to in- 
crease \Mnr\ to 4.3 would result in i? 42, just above the upper limit, but the resulting V 



would remain compatible with the bounds in Eq. (1.1). 



Attempting instead to decrease |M^^| below 4 results in a mass-squared difference ratio that is 
too small. For example, |M^t-| = 3.9 would result in R ^ 23, which is less than the empirical 
lower bound of 25.5. As for the previous case, the resulting V would be compatible with data. 



Thus in the //r-symmetric case, the experimental constraint on R 



m 



31 



/m2i is more strin- 



gent than the constraints for the entries in V. This can be understood from Fig. [T| in which 



we plot i? as a function of |M^t-|/M^^ for the matrix in Eq. (5.2). As |M^t-|/M^^ — )• 1 we 



have m^i 



— )■ 0, so that the ratio R diverges as \M. 



detail in the next section (see Eq. (6.1)). 



I1T\ 



— )• M^^, which we will discuss in more 
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10 : 

OjO 0.4 0.6 OA ]jj K^^TlZ-fJ/J-^ 

Figure 1: The mass-squared difference ratio R = rn^i/m^i for the ^r-symmetric 
ansatz Eq. ( |5.2[ ). The curves correspond to the fixed values M^^ = 3 (Red), 4 
(Green), 5 (Gray), 6 (Blue), 7 (Black). For |M^^|/M^^ > 1, reflect the graph about 
the vertical line |M^t-|/M^^ = 1 




\M MM 



In comparison, there is no divergence in either Vei or V^2, whose sensitivity to the ratio 
\Mfj_r\/ Mfj^^ is displayed in Figs. 2(a) and 2(b) 



iVell 

0.95 
0.90 
0.85 
0.80 

0.75 

0.™ 

0.65 1 

0.0 





(a) iKil vs. |AV|/M^^ (b) vs. |7VV|/Af^^ 

Figure 2: The elements Vei and Ve2 of the mixing matrix for the /xr-symmetric ansatz Eq. 
(5.2). The curves correspond to the fixed values M^^ = 3 (Red), 4 (Green), 5 (Gray), 6 (Blue), 



7 (Black). For \M^r\/M^^ > 1, reflect each graph about the vertical line |M^t-|/M^^ = 1. 



Having gained an analytic understanding of the mass matrices in Class I (3.7) and Class II 
(3.8), as well as their yur-symmetric intersection (5.2), we now turn to numerics. The analy- 
sis that follows should be useful for classifying perturbations away from tribimaximal mixing 
within the /ir-symmetric ansatz as well as for classifying deviations from /ir symmetry in 
more general mass matrices. 
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6 Real Mass Matrices: Class I {Me^, = Mer) 



We now begin a numerical study of the Class I ansatz defined by Eq. (3.7), which for the 
convenience of the reader we display again: 



^0 1 



1 



Mi = I M^^ 

Mr 

Here Mnr is strictly negative, and all other nonzero entries are strictly positive. 



Figure 3(a) shows a plot of the allowed values for the ratio M^^/M^^ while letting M^ 



range over all its possible values. We find that the nonzero diagonal entries can lie in the 
ranges M^^ ~ 2 — 9 and Mrr ~ 2 — 10. The fact that these ranges are essentially the same 
is something we already knew, since as discussed in Section [3] the case M^^ = Mrr is the 
/xr-symmetric subcase of Class I. 



Mr 



(a) M 




T-T- VS. M^^J^ for all allowed values of |Af^T-| 




(b) 



\M^r \ vs. M^^ with Mrr 



M„ 



Figure 3: The ratio M^^^/Mrr in Class I {Me^j. = Mer = !)• Figure 1(a) shows the values of 
the diagonal entries in Class I for all possible allowed values of M^r- Figure 1(b) shows the 
allowed values for M^r for matrices of Class 1 that also satisfy M^^ = Mrr and therefore are /ur- 
symmetric. As discussed in Section[5| mass matrices with M^^ = Mrr for which |M^^| = M„„. — 1 



(lower solid line) and = M^,- + 1 (upper solid line) yield tribimaximal mixing. 



^"Here and throughout the rest of the paper, we use the "~" symbol to denote a rough guide for the values 
of the entries in M, to be compared with either = M^r = 1 (Class I) or = Mrr = 5 (Class II). 
The idea is to get a feel for what the entries in M can be, and then afterwards to hunt for precise numerical 
values that fit data. 
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However, looking at Figure [3(a) in isolation may give the misleading impression that the case 
Mt-t is allowed when in fact it is experimentally ruled out, as can be seen 



numerically in Figure 3(b) This can also be seen in Figs. [T| |2(a) and |2(b) when compared 
with the bounds given in Eqs. (1.1) and (1.2). 



This can also be understood analytically as follows. The mass matrix 

^^0 1 1 

Mb 



a —a 
a 



(6.1) 



implies a "bimaximal" mixing matrix 



PMNS — 



1 1 

2 2 



11-1 



(6.2) 



\ 2 2 v^/ 

and two equal neutrino masses, both of which are incompatible with the empirically allowed 



ranges quoted in (1.1) and (1.2). 



Figure 3(b)| shows that the set of allowed mass matrices splits into two branches, with 
larger and smaller |M^t-|, which yield a larger and smaller R respectively. For example, for 

5 we recoveirn either the matrix with M,,^ = —4 and i? 27 in Eq. 



fixed M, 



(3.3), or the matrix with M. 



-6 and i? = 40 in Eq. (3.5). 



We can use Figures 3(a) and 3(b) to look for examples of mass matrices with non-tribimaximal 



mixing. Towards the upper limit of |M^^| ~ 8, we find 

M= I 6 -7.6 1 ^ R^35 and |V"pmns| 




0.86 0.50 
0.36 0.61 0.71 
0.36 0.61 0.71 



Like tribimaximal mixing, this case has 



Me, and M^^ 



M^^ with Vr 



a3 



(0, 



-1 



Unlike tribimaximal mixing, the second column of V is not proportional to (1, 1, 1), and the 
first column changes accordingly to maintain orthogonality. This is all consistent with the 
analytic understanding of the /ir-symmetric ansatz from Section [5j 

Thus, as emphasized throughout, this is an example with M^^ = Mrr but without tribi- 
maximal mixing. In passing, we mention that increasing |M„,| to 7.7 would make \Vei\ too 



large and Ve2 too small with respect to the bounds given in Eq. (1.1). The sensitivity of V^i 
and Ve2 to changes in |M^t-| can be seen in Figs. 2(a) and 2(b) 



Toward the lower limit of |M^,| ~ 3, we find 

'0 1 1 ' 
M = I 3 -2.9 
3.9 



i?^27 and IK 



PMNS 



0.77 0.64 0.02 
0.48 0.59 0.65 
0.42 0.49 0.76 



As discussed in Section 3l we also have Eq. (3.5) 
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which is of a similar form to the "nTB" matrix of Eq. (3.6), except with M-rr larger than 
Again in passing, we point out that increasing |Mur| to 3 would make \Vei\ 



M, 



MM 



\M, 



too small and Ve2 too large. On the other hand, keeping |M^t-| = 2.9 but increasing Mrr 
to 4 would result in R 25.2, just below the experimental lower bound, while maintaining 
a consistent mixing matrix V . An example of the sensitivity of R to the ratios of various 
entries in M can be seen in Fig. [TJ although the matrix above is not yur symmetric. 



Numerically we find that for all allowed values for M^^ and Mrr-, the entry 



\M^r\ can be 



in the range ~ 3 — 8. This can be seen in Figs. 4(a) and 4(b) 



\M.r\ 



\M.r\ 




M 




Mr 



(a) |Af^tT-| vs. M^^ for all allowed AIt-t (b) iM^T-j vs. Mrr for all allowed A/^^ 
Figure 4: The allowed values of \M^r\ in Class I {Me^ = Mer = 1). 



We conclude the study of Class I with a comment on VJ=3. Figures 5(a) and 5(b) show that 
these mass matrices exhibit a maximum value |Ve3| ~ 0.03, which is rather small. However, 
these figures also identify that having M^^ and Mrr both less than 3 or greater than 8 ensures 
a nonzero Ves. Recall that previously we observed M^^j ~ 2 — 9 and Mrr ~ 2 — 10, so that the 
narrow ranges M„„ ~2 — 3or8 — 9, and Mrr ~ 2 — 3 or 8 — 10 are those which necessarily 



produce V^s ^ 0. These ranges are correlated, so that M^^ 
As an example, we have: 



2 with Mr 



10 is not allowed. 







M 




R^29 and Vt 



PMNS 



and thereby generate |V^3| ~ 0.02, as for MnXB- 



-0.86 0.51 0.02 
0.43 0.70 0.58 
0.28 0.50 -0.82 



In summary, matrices of Class I necessarily have "small" values of V^s, reaching a maxi- 
mum of only ~ 0.03. 



12 





(a) Ve3 vs. Afpp for all allowed Mtt and |A^^t| (b) Ves vs. Mrr for all allowed Af^^, and |Af^,-| 



Figure 5: VeS as a function of Afa^ in Class I (Afg/x = M^r = 1). In both plots, the 
variables not displayed explicitly on the axes are allowed to range over all of their possible 
values that result in an acceptable mixing matrix V and mass-squared difference ratio R. 



7 Real Mass Matrices: Class II (Af^^ = Mrr) 

Now consider matrices from Class II, which is defined by 



'0 Mef, M,, 



Mil 




In Figures 6(a) and 6(b) we plot the allowed values of IM^.^! as a function of Mg^ and 



MgT-. We find that |M^t-| is constrained to be very close to either 4 or 6, reminiscent of the 



tribimaximal cases given in Eqs. (3.3) and (3.5). This tells us that matrices in Class II 



necessarily exhibit near tribimaximal mixing if the ratio Mg^/Mer is close to 1. 



IM„,| 



\M.r\ 



itttMHIMIimilllMIt* 



M, 



(a) |AfpT-| vs. for all allowed values of Aler (b) |Af^T-| vs. Afer for all allowed values of Afg^ 

Figure 6: The allowed values of \Mfj,r\ as a function of Afg^ and Afer for Class II (Af^^ = M^r = 
5). Parameters not displayed explicitly on the axes are allowed to attain all values compatible 



with (1.1) and (1.2). 



We emphasize that the reason the ratio Mg^/Mer characterizes proximity to tribimaximal 
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mixing in Class II is simply because the data in (1.1) and (1.2) constrain |M^^ 
to 4 or 6 (in units for which M^^ = Mrr = 5). Otherwise, as discussed below Eq 



to be close 



(5.2) 



symmetry does not imply tribimaximal mixing. 



In Figure It] we examine the ratio Mi.^/Mer- The case Mg^/Me,- = 1 is allowed for M, 



Mer ~ 



0.8 — 1.5, but for values out of this range for either Mg^ or M^r, the mixing matrix 
will deviate significantly from the tribimaximal ansatz while still fitting data. 



<•••••• 

<••••••• 



M. 



e/i 



Figure 7: The ratio Mer/Me^, (for all allowed \M^r\) in Class II (M^^ = Mrr = 5). 



In particular, it is possible for matrices of Class II (M^ 



with either Men = or Mer = but not both. For example: 



M 








i?^41 and |V^pmns| 



Mrr, ^ M, 



0.50 0.11 
0.66 0.69 
0.56 0.71 




to fit data 



Decreasing Mer to 1.3 would make R and |Vei| too large and Ve2 too small with respect to 
the bounds in Eq. (1.1 ). This two-zero texture is labeled "Case Ai" in a study by Frampton, 



Glashow and Marfatia about possible zeros in the neutrino mass matrix in the flavor basis 
[T8] . The salient feature of both their work and ours is the possibility of a "large" Vgs, on 
which we now elaborate. 



For either Mg^ or M^ 



less than ~ 0.6 or greater than 



1.6 in Class II, a nonzero 
is generated. Figures 8(a) and 8(b) show that the empirical upper limit \Vez\ ~ 0.2 can be 
generated for Mg^ or Mer close to ~ or 3. Using these along with Figure [7], we find: 







M 




R^33 and |Vf 



PMNSI 



0.79 0.57 0.22 
0.32 0.69 0.65 
0.52 0.44 0.73 



Decreasing Mg^ to zero would generate | Ksl ~ 0.23, which is too large if we beheve the upper 



bound given in Eq. (1.1). 
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Ve3 



Ve3 



Si0 



i 





M. 



(a) Ve3 vs. Afe^ for all allowed AI^t and |M^^| (b) Ves vs. Met- for all allowed Mg^^ and |A/^t-| 
Figure 8: Ves as a function of Mga in Class II (M^^ = M^,- = 5). 



In summary, mass matrices of Class II can result in values for |Ve3| anywhere from to 
the empirical upper limit of ~ 0.22. In particular, when the ratio Me^/Me,- is greater than 
~ 2 or less than ~ 1/2, a "large" jVesI is necessarily generated. 

8 CP Violation 

We will now allow for the possibility that the neutrino mass matrix violates CP. As discussed 
in Section |4| any 3-by-3 complex symmetric mass matrix with M^e = can be rephased into 



the form of Eq. (4.2), which we repeat for convenience: 



'0 |MeJ \M,r\e'^' 
M= I \M„„\ \M,Je'^ 



TT 



We need to diagonalize M to determine how cp and rj contribute to the CP-violating angle 



Sep and to the Majorana phase angles p and a. (Recall the notation of Eq. (1.3).) The 
mapping of two phases (p and r] to three observables 6cp, p and a is explained by the fact 
that p and a are not independent parameters when Mee = 0. With a complex mass matrix, 
the condition Mee = implies 

tan^ 9, = -( ole-''^^ + 4e-^-mA 



which is the generalization of Eq. (5.1) with the possibility of p and a being different from 



or tt/2. The imaginary part of this fixes p in terms of a through the relation 

^ = -tan^^3^ (8.2) 
sin(2cr) mi 

so that only one of these phases is an independent parameter. 

The main result of the generalization to complex mass matrices is that nontrivial phases 
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open up new regions for the allowed values of |M, 



To understand this claim it is sufficient to specialize to the following example: Recall the 



matrix Mb from Eq. (6.1), which predicted mi = m2 and bimaximal mixing with |Vei| = Ve2, 



which are incompatible with the bounds in (1.1) and (1.2). Generalizing this matrix to the 
complex case 

^0 1 e' 



Ml 



Be 



a ae 
a 



IT] 



1.3) 



can split the degeneracy mi = 1712 and modify |V^i| = Ve2 significantly enough to become 
compatible with oscillation data. In the next section, we will show that the matrix Mbc 
can result in tribimaximal mixing with rrii = (1,2,9) and rrii = (1,2,11), just as in the 
CP-conserving case discussed in Section |3} 



9 Tribimaximal Mixing and Nonzero Majorana Phases 



Consider tribimaximal mixing^, meaning H'mns = Vtb, but with arbitrary Majorana phases 
so that^^ M = V*DV'^ is complex even though Sep drops out since V^3 = 0. 



For tribimaximal mixing with arbitrary Majorana phases, the condition Mee = now fixes 

2mi e"^^" + ma = (9.1) 



which corresponds to taking 62 = and 6^ = sin ^(l/y/S) in Eq. (8.1). This implies 

1712 = 2mi and p = a + [2n — 1) — 



(9.2) 



where n is any integer. Consider the case p = —a = — 7r/4. Upon choosing the Ki as given 
above Eq. (4.2), we find a rephased mass matrix M of the form 



/O 1 



M = Mtbc 



a ae 



vq 



(9.3) 



where 



+ ^ and = TT + tan ^ (^-|^^ . 



(9.4) 



As displayed above, this matrix has ^9 = 0. We have set the overall scale rriy = nii = M^fj, 

Mer to 1. 



^^Some of our work in this section overlaps with that of Z. Z. Xing f9]. 

i3We remind the reader that D — diag(mi, m2, 771,3) is real and positive, and V = /CVpmnsA^ is the full 



3-by-3 unitary matrix including the extra phases in K. and Ai. If Vi denotes the 



■th 



'PMNS-/ 

column of V, 



then 



^■^As discussed below Eq. (4.3), we can equivalently set ip — tt ii we subtract tt from rj. 
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For (mi, 777,2) = (1,2) we can invert the definition R = in\i/m\i to get 777,3 = v^3-R + 1 and 
thus obtain the parameters a and 7/ as a function purely of the experimentally constrained 
ratio R: 

1 - ^ 7/ = tan-i i^^^R+A (9-5) 



V3^ 



and 



where 25.5 < R < 41.0 as given in Eq. (1.2). These parameters are plotted as a func- 
tion of R in Figs. 9(a) and 9(b) In units of mi = 1, the mass of the heaviest neutrino is 
8.80 < ms < 11.1. 



a 



(a) a vs. R for tribimaximal mixing 



R 



TJ-TT 



R 



(b) rj — IT vs. R for tribimaximal mixing 



Figure 9: The parameters a = 



Mr 



Lfj_r\ and rj = ai'g(Mf^r) in the matrix Mbc given 
in Eq. (8.3), for the particular case in which the mixing matrix is exactly tribimaximal. In this 
case we have m2 = 2mi and Mg^ 
7772 = 2 and thus 8.80 < 7713 



ier — "^1) SO that setting Mg^ 
VSR + 1 < 11.1. 



Mer = 1 implies 7711 = 1, 



Consider the two examples ms = 9 and ms = 11. For m^ = (1,2,9) we have a ^ 4.5 and 
Tj — TT K, 0.22, and for m^ = (1, 2, 11) we have a ~ 5.5 and rj — n ^ 0.18. This interpolates 
between the two branches |M^t-| ~ 4 and |M^t-| ~ 6 of the real- valued tribimaximal mass 
matrix. For instance, m^ = (1, 2, 10) (so that R = 33 exactly) implies a = | vToT ^ 5.0 and 
r]-TT ^0.20. 

Therefore in addition to the CP-conserving case with (p = 0, rj = it and \M^r\ 7^ M^fj^, 
we find a new class of allowed tribimaximal mass matrices with (f = 0, rj ~ 7r + 0.2 and 
\M^T-\ = M^fj^. It is important to note that, in contrast, the case rj ~ 0.2 is not allowed 



unless the angle cp is changed to ~ tt. (Recall the notation of Eq. (8.3).) This should be 



understood in the context of the discussion below Eq. (4.3), in which we showed that phases 



of e^*'^ = — 1 can be exchanged between Me^ and M„ 



10 Complex Mass Matrices with M^^ = iMer] and M^^ = Mrr = \Mf,r\ 



We now turn to a numerical study of the matrix Mbc given in Eq. (8.3). Up to the phases, 
this is the /ir-symmetric subcase of both Classes I and II with the additional condition 
\M^T-\ = ^nn- (Note that for non-tribimaximal mixing, the phase (p is no longer necessarily 
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zero or 7r.) An immediate striking feature of this matrix is given in Figure |10(b)| which 
shows that Ir^l can take essentially only two possible values: ~ 0.2 and ~ n. This corroborates 



the intuition we gained from tribimaximal mixing with |M„t-| = M. 




(a) tp vs. for all allowed values of rj (b) rj vs. M^^ for all allowed values of (p 



Figure 10: The phase angles as a function of M^^ 
given in Eq. 



Marl the matrix Mbc 



3). The values near r] 



±0.2 in (b) should be understood in the context 
of the discussion below Eq. (4.3). That is, the allowed values {ip,r]) = (0,±(7r + 0.2)) 
are phenomenologically equivalent to the values (v^, ??) = (vr, ±0.2). In contrast, the values 
{ip,r]) = (0, ±0.2) are not compatible with the data in Eqs. (1.1) and (1.2). 



Figures 11(a) and |ll(b) show that this CP- violating matrix interpolates between the real- 



valued cases with |M^^| = 4 and |M„^| = 6, which give R ^ 27 and i? = 40 respectively. 



R 





R 



(a) R vs. (f for all allowed 77 with Af^ 



Figure 11: The values of R for the particular cases M^^ 

6 in the matrix Mbc of Eq, 



(b) R vs. ip for all allowed 77 with M^^ — Q 



|M, 



3 and 



3). The angle t] ranges over all 

0.2 



allowed values, which as shown in Fig. |10(b)| amounts to only the possibilities r] 
(with ~ vr) and 77 ~ vr (with ~ 0). 



Recall that the matrix Mbc of Eq. (8.3) for the special case (p = and rj = n reduces to 
the matrix Mb of Eq. (6.1), which implies a bimaximal mixing matrix and thus |V"ei| = Ve2, 



which is incompatible with the bounds given in (1.1). Figures 12(a) and 12(b) show 
that complex phases can generate mixing matrices that fall in the empirically allowed range 
0.14 < iKil - Ve2 < 0.36. 
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|Vel|-ye2 



|Vel|-ye2 



(a) |V"ei| — Ve2 VS. (fi for all allowed rj with A/^^ — 3 (b) \Vei \ — Ve2 vs. (p for all allowed 77 with M^^ — 6 



'MM 



\M, 



tJ.T\ 



3 and M, 



MM 



Figure 12: The values of |T4i| — V^e2 for the particular cases 
Mrr = l-^^rl = 6 iu the matrix Mbc of Eq. (8.3). Recall from Eqs. (6.1) and (6.2) that in the absence 
of complex phases in Mbc, the mixing matrix is of bimaximal form with \Vei\ — Ve2 = 0, which is 



experimentally ruled out. The angle r/ ranges over all allowed values, which as shown in Fig. 10(b) 
amounts to the two possibilities r/ ~ 0.2 and ~ vr. 



X 



|Ve3| 



y 



_i I L_ 




(a) Ve3 vs. (/3 for all allowed 77 with M^^ = 3 
Figure 13: The value of Ves for the particular cases M^^ 
6 in the matrix Mbc of Eq. (8.3). 



(b) Ve3 vs. Lp for all allowed t] with Af^^ = 6 



M^r = WaA = 3 and M, 



M. 



MM 



The angular parameterization of the mixing matrix makes clear that 5cv only contributes to 



neutrino oscillations when Vez 7^ 0. Figures 13(a) and 13(b) show that the magnitude of 
Vez depends strongly on the value of a = \M,,A = M, 



MM 



M^-T- in mass matrices with Mp, 







and M 



iMeJ = 1. 



More generally, amplitudes for CP-violating oscillation processes are proportional to the 
rephasing-invariant quantity J = — Im(\43V^2V'g2^3) [IS]- Figures 14(a) and 14(b) show 
that, like V"e3, the quantity J also depends strongly on the value of a. For a near its lower 
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bound of ~ 3, the quantity J is of order ~ 10^^, but for larger a ~ 6 we find that J is at 
most ~ 10^^ and can drop to zero. For an example with a large jVesI and nonzero 5cp; and 
hence a large J, we find: 



M 




R^39 and IK 



PMNSI 



0.79 0.58 0.19 
0.52 0.50 0.70 
0.32 0.65 0.69 



with 6cp 
we find 



—0.38 and J !^ 1.6 x 10 ^. In contrast, by increasing |M^,-| = M^^ = M^-t to 6 



M 




R^39 and IK 



PMNSI 



0.83 0.55 0.02 
0.40 0.58 0.71 
0.38 0.60 0.71 



with 6cp ~ —2.7 X 10~^ and J ~ 1.3 x 10~^. These two examples were chosen intentionally 
to yield the same value for R. Note that for \M^r\ = = Mrr = 6, keeping ip = 1.1 
would result in i? ^ 64 with \Vei\ ~ 0.90, both of which are too large, and in Ve2 ~ 0.43, 
which is too small. On the other hand, the other entries in Vpmns would all stay within the 
empirically allowed ranges. 

In summary, we learn that for the matrix Mbc decreasing the value of M^^j,/Me^ increases 
the value of iKsl and thereby results in the possibility for larger amplitudes for CP- violating 
processes. 



_i I I i_ 



V 



_l I I I I l_ 



(a) J vs. f for all allowed 77 with Af^ 



J 



* 

* 



if 



Figure 14: The value of J 
3 and M^^ = Mrr = \M^r\ 



(b) J vs. (/5 for all allowed 77 with M^^ = 6 



lvti{Vei,V^2V*2V*.^) for the particular cases M^^ = Mrr = 



6 in the matrix Mbc of Eq. (8.3). 



^1T\ 
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11 Discussion 



We have studied three types of neutrino mass matrices in the flavor basis with M^e = 0. 
The first two types are the CP-conserving matrices of Class I (Mg^ = M^.^.) and of Class II 



(M^^j = M^T-), which we display again for the convenience of the reader (see (3.7) and (3.8)): 



'0 



Ml 



M„„, M, 



tJ.T 



and Mil 



'0 Me, 



fll.ll 



The intersection of these two classes is the /ir-symmetric ansatz (see (5.2)) 



fiT-sym 




(11.2) 



which as discussed should be thought of as a useful phenomenological starting point. 



The salient phenomenological distinction between Classes I and II is that mass matrices 
of Class I can accommodate only a small V^3 up to ~ 0.03 (Figs. 5(a) and 5(b)), while mass 
matrices of Class II can predict an arbitrarily large Ves (Figs. 8(a) and |8(b)' ). Thus fundamen- 
tal theories which predict a neutrino mass matrix with M^^ ~ 



Me^ > Mer 

of Ks- 



as opposed to Men ~ M, 



Mrr and either Me^ -C Mer or 
will be the most constrained by future measurements 



The third type of matrix we studied is the complex matrix (see (8.3)) 

^0 1 e^^ 



Mbc = 



a ae 



(11.3) 



For this matrix, smaller values of a result in larger values of and J = — Im(V^3V^2K2V^*3), 
and thus provide experimentally promising signals of CP violation in neutrino oscillations 



(Figs. 13(a) and 14(a)). In contrast, larger values of a drive Ves and J to zero (Figs. 13(b) 
and|14(b)D. 



A particularly interesting example is obtained from Mbc for the particular case = with 



1 



a=-\/m?^ + l and 77 = vr -|- tan ^. „ 

2 V ^ \ml-l 



21713 



(11.4) 



Here 777,3 is the mass of the heaviest neutrino in units of the lighest neutrino (777,3 > ^2 > 1^1), 
and the other two masses are (7771,7772) = (1,2). In this case the mixing matrix is exactly 
tribimaximal, even though M^^ = M^t = \Mfj_r\. (See Section |9 
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